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We present X-ray observations of the as of yet unidentified very high-energy 
(VHE) 7-ray source HESS J1640-465 with the aim of establishing a counterpart 
of this source in the keV energy range, and identifying the mechanism responsible 
for the VHE emission. The 21.8 ksec XMM-Newton observation of HESS J1640- 
465 in September 2005 represents a significant improvement in sensitivity and 
angular resolution over previous ASCA studies in this region. These new data 
show a hard-spectrum X-ray emitting object at the centroid of the H.E.S.S. 
source, within the shell of the radio Supernova Remnant (SNR) G338. 3-0.0. 
This object is consistent with the position and flux previously measured by both 
ASCA and Swift-XRT but is now shown to be significantly extended. We argue 
that this object is very likely the counterpart to HESS J1640-465 and that both 
objects may represent the Pulsar Wind Nebula of an as of yet undiscovered pulsar 
associated with G338. 3-0.0. 
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Introduction 



The last years brought an emerging new population of Galactic VHE 7-ray sources. More 
than 20 new objects are now established as sources of 7-rays above 100 GeV by Atmospheric 
Cherenkov Telescopes such as H.E.S.S. (jHintonll2004l ). In particular a first unbiased survey 



of the Inner Galaxy by H.E.S.S. found more than a dozen new sources ( jAharonian et al. 



2005al . l2006al ). While several Galactic 7-ray s ources can be identified with counterparts 
at o ther wavebands, such as RX J1713.7-3946 (lAharonian et all booi l2006bh or MSH-15- 
52 ( lAharonian et al.l l2005bT). most of them have not yet been identified with a striking 
positional counterpart (jAharonian et al.ll2006al ; lFunkll2006l ). A counterpart search proceeds 
in several steps: first, a positional counterpart at other wavebands must be established, typ- 
ically through high-angular resolution observations in the radio and X-ray bands. In these 
bands non-thermal particle populations can be studied through their synchrotron emission 
without being strongly absorbed by interstellar dust. Following the detection of an astro- 
metric counterpart, a firm identification requires both a viable 7-ray emission mechanism 
and a consistent multi-frequency picture of the positional counterpart. Here we report on 
X-ray observations on HESS J1640-465 with the X-ray XMM-Newton satellite as part of an 
ongoing programme to study unidentified VHE 7-ray sources with X-ray satellites. 



H ESS J1 640-465 was first detected by H.E.S.S. in a survey of the inner Galaxy (jAharonian et al 



2005al . l2006al ) in 2004. Following the initial discovery the source was re-observed in pointed 
observations with a resulting statistical significance of ~ 14a in 14.3 hours of lifetime, with 
a total of 313 7-ray events above the background of ~ 300 events. The 7-ray emission profile 
is Gaussian and rather compact in nature with an rms width of 2.7' ± 0.5', less extended 
than most of the sources detected in the survey. The 7-ray energy spectrum was reported to 
follow a simple power- law with photon index 2.42 ± 0.15 and a total integrated flux above 
200 GeV of (20.9 ± 2.2 stat ) x 10~ 12 cm" 2 s" 1 (corresponding to 2.2 x 10~ n erg cm" 2 s" 1 ). 
HESS J1640-465 is located in the Galactic plane at an RA-Dec 20 oo of 16 h 40 m 44 s , -46 d 31'57" 
(or I = 338.32°, b = —0.02°). This direction on the sky lines up with the prominent 3-kpc- 
arm-tangent region of the Galaxy at a distance of 8 kpc. 

HESS J1640-465 was found to sho w compelling positional coincidence with the shell-type 
Supernova remnant (SNR) G338. 3-0.0 (jGreenll2004j ). This broken-shell SNR with a diameter 
of 8' was detected i n the 843 MHz rad io survey using the Molonglo Observatory Synthesis 
Telescope (MOST) (jGreen et al.lll999l ). G338. 3-0.0 is located at a RA-Dec 2 n nn position o f 
16 h 41 m 00 s , — 46 d 34' and is not particularly well studied by radio telescopes (jGreenl 12004 ). 
A reanalysis of the MOST data (T. Cheung, private communication) yields a total radio 
emission from the remnant of ~ 10 Jy on a background of ~ 10 mJy per beam outside of 
the remnant. The central point source within the remnant shows a flux of ~ 160 mJy with 
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a peak flux of ~ 100 mJy. The connection between this central source and the remnant is 
as of yet not evident. The radio emission from the interior of the remnant away from the 
central point-source amounts to ^20 mJy. No pulsar i s known within G338. 3-0.0 in the 



current ATNF pulsar catalogue ( Manchester et al. 2005h. The proj ection of G338. 3-0.0 is 



located at the edge of a bright Hll-region (IWhiteoak &: Greenlll996l ) (located at a distance 



of 3 kpc) with ongoing star formation and dust emission. This region seemingly conne cts the 
radio emission of G338. 3-0.0 with the close-by Green's catalogue SNR G 338.5+0.1 (jGreen 



20041 ). The best distance estimate of G338. 3-0.0 so far comes from the relation of surface 



brightness to distance which places G338. 3-0.0 at a distance of 8.6 kpc, right within the 
3-kpc-arm-tangent region (and leaves it thus unrelated to the HII- region). At X-ray energies 
the region has been observed in the ASCA Galactic plane surve y, resulting in the de tection 
of AX J1640. 7-4632 in positional coincidence with the SNR (ISugizaki et al.l l200ll ). The 
ASCA source showed a rather soft photon index of 3.0l ;g, a total unabsorbed flux between 



2 keV and 10 keV of F 2 - W ke v 



1.2 x 10 12 erg cm 2 s 1 , 



and a rather high absorption of 



nn = 9.6^3 3 x 10 22 cm~ 2 com pared to the total column density in this direction (ncaiH = 
2.2 x 10 22 cm -2 as reported by iDickey fc Lockmanl (119901 ) ). Recent Swift XRT observations 
on HESS J1640-465 revealed one X-r ay point-source wi thin the radio shell and in coincidence 
with AX J1640. 7-4632 (source SI in lLandi et ap200fih . The XRT data confirm the spectral 
parameters found by ASCA resulting in still rather unconstrained values of the photon index 
T = 2.6+i;? and the 2-10 keV flux F 2 - 
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0.7 x 10 12 erg cm 2 s l . 



At 7-ray energies above 100 MeV, 3EG J1639-4702, an unidentified EGRET source ( IHartman et al. 



19991 ). is found spatially compatible with HESS J1640-465 at an angular distance of 34'. 



3EG J1639-4702 was only accepted for inclusion in the 3EG catalogue based on its cumulative 
significance over 4 years of EGRET observation, but never satisfied the detection threshold 
in any individual EGRET observation. Consequently, there is no reliable variability estimate 
available for this source, and a classification based on the available EGRET observables ap- 
peared impossible. 3EG J1639-4702 was labelled as an extended and confused 7-ray source, 
which puts the rather tight positional coincidence of 34' with HESS J1640-465 somewhat 
into perspective. Various positional counterparts were found: a rather specula tive suggestion 



was t o associate it both with the unidentified ASCA so urce AX J1639.0-46 42 (ISugizaki et al. 



200l|), and the INTEGRAL source IGR J16393-4643 (IMalizia et al.l 120041) conc luding that 



3EG J1639-4702 is a dust-enshrouded microquasar candidate ( jCombi et al.ll2004l ). This sug- 
gestion was revised by the detection of a heavily obscu red X-ray pulsar in XMM-Newton and 



RXT E observations coincident with IGR J16393-4643 (IBodaghee et all2006l ; iThompson et al. 



20061 ). The corresponding mass estimates allowed for a classification of IGR J16393-4643 as 



a high-mass X-ray binary system. AX J1640. 7-4632 and IGR J16393-4643 might therefore 
not be related to 3EG J1639-4702 given that this would constitute the first association be- 
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tween an unidentified EGRET source and a highly obscured INTEGRAL-detected HMXB. 
This scenario currently lacks both viable multi-waveband modelling and predictions for GeV 
energies. Other associations for 3EG J1639-4702 were suggested, including posi tionally coin- 
cident radio pulsars (PSR J1637-4642, PSR J1640-4648, and P SR J1637-4721) JTorres et al. 



20011) and SNRs (G337.8-0.1, G338. 1+0.4, and G338.3-0. 0l dRomero et al. 



199j3). Two of 



these pulsars can be disregarded on energetics grounds (jTorres et al.l 120031 ). leaving only 
PSR J1637-4721 with a moderate requi red efficiency of 12% for conversion of spin down 
power into EGRET 7-radiation. Finally, iRomero et al.l (119991 ) also listed the stellar cluster 
and OB-associations NGC 6204, and Ara 1A A as positional coincidences, the latter addition- 
ally coinciding with 3EG J1655-4554. Summarising the possible associations of 3EG J1639- 
4702, the typical situation of an observationally weakly constrained EGRET detection allows 
for a variety of interpretations. A unique counterpart association will only be possible if the 
source location uncertainty is refined (using GLAST-LAT), or intriguing features at other 
wavebands are discovered, as perhaps in this case from H.E.S.S. observations. 

The positional coincidences of HESS J1640-465 with G338.3-0.0 and AX J1640. 7-4632 
(and possibly also 3EG J1639-4702) suggests a connection between these objects. However, 
from the described multi-frequency data it remains so far unclear whether the X-ray and 
7-ray emission originates within the shell of the radio SNR or rather in connection with 
a central X-ray object embedded within the shell. If the distance estimate of 8.6 kpc is 
correct, the radius of the SNR is ~ 10 pc and the Sedov solution for the evolutionary stage 
of the SNR yields ages between ~ (0.2 — 2) x 10 4 years (assuming a typical ambient density 
between ~ 0.1 — 10 cm -3 and an explosion energy of E = 10 51 erg). In particular for 
densities larger than the canonical lcm -3 , the Sedov solution yields an age in excess of the 
typical VHE-bright 7-ray SNRs, such as RX J 1713. 7-3946 and RX J0852.0-4622 with ages 
of ~ 2000 years or younger. In this letter we report on X-ray observations towards the 
G338.3-0.0/HESS J1640-465 region performed with XMM-Newton with the aim to study 
the connection between the radio/X-ray and 7-ray emitting regions in the light of higher 
angular resolution X-ray data. 



2. XMM-Newton Observations of the region 

XMM-Newton observed HESS J 1640-465 on the 20th of August 2005 for 21.8 ksec in 
satellite revolution 1043 (observation id 0302560201). All the instrument cameras (EPIC 
MOS1, MOS2, PN) were operated in full-frame mode with a medium filter to screen out 
optical and UV light. For data analysis, the XMM-Newton Science analysis software (SAS) 
version 7.0 was used together with the Extended Source Analysis Software package (XMM- 
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ESAS) version 1.0 (ISnowden. Collier fc Kuntzl l2004j ). Following standard data reduction 



and calibration procedures, the data set was cleaned from temporally occurring background 
caused by soft proton flares. The resulting observation time amounts to only 7.3 ksec of 
useful data. Figure [T] shows an adaptively smoothed count map below 2 keV (top) and 
above 2 keV (bottom) combining events for the MOS-1 and MOS-2 detectors of XMM- 
Newton of the region surrounding HESS J1640-465. The white contours show the 843 MHz 
MOST-detected radio emission. Three prominent X-ray sources are visible, coinciding with 
the recently reported Swift XRT sources (and accordingly labelled #1, #2, and #3). Source 
#1 positionally coincides with the ASCA source AX J1640. 7-4632 and is located within the 
radio shell of G338. 3-0.0. No shell-like X-ray emission such as seen in the MOST data is 
apparent in the XMM- Newton data set. Also visible is the band of stray-light crossing the 
field of view from the north to the east, which can be attributed to the close-by (~ 0.7°) low 
mass X-ray binary GX 340.0. This stray-light contribution to the spectral analysis within a 
small region surrounding the X-ray source AX J1640. 7-4632 seems to be rather small, since 
various background estimation methods using regions within the same field of view or from 
blank sky observations yield results comparable within the statistical uncertainties. 

In a first step to characterise the X-ray sources in the region, the standard source 
detection algorithm emldetect has been applied to the data in the energy range from 0.5- 
10 keV, and as a cross-check in the restricted energy bands 0.5-2 keV, 2-4.5 keV, and 4.5- 
10 keV. Three X-ray sources are detected as depicted in Figure [1] by the cyan circles. The 
detection of source #1 in the 0.5-2 keV energy band is very weak, indicating a rather hard 
source with a high degree of absorption at low energies. Source #3 on the other hand is only 
detected in the 0.5-2 keV band, indicating a rather soft source. While these t hree sources 



are lo cated within the several-arcsecond error bars of the Swift XRT-sources ( jLandi et al. 



20061 ). the faint XRT-sources #4 and #5 are not detected in this XMM-Newton data, neither 
in the MOS-detectors in which they fall into gaps between the CCD chips nor in the PN- 
detector in which they are covered by the stray-light from GX 340.0. Table [T] summarises the 
properties of the three sources for the 0.5-10 keV energy band. Only source #1, coincident 
with the ASCA source AX J1640. 7-4632 (in the following named XMMU J164045.4-463131) 
is extended in nature, consistently in all energy bands. Figure [2] shows a zoom on the region 
surrounding XMMU J164045. 4-463131 along with the best fit po sitions for the XMM-Newton 



data presented here and the Swift XRT data (jLandi et al.ll2006l ). The source position fitting 



tool emldetect determines the extension of XMMU J164045. 4-463131 to be incompatible 
with a point-source. Using a Gaussian model, the extension in the best-constrained energy 
band above 4.5 keV is determined to be 10.7 ±0.1 image pixels, corresponding to a width 
on the sky of 27" ± 3". For the lower energy bands, the best-fit width is larger but may be 
affected by stray-light, the band 4.5 keV to 10 keV was therefore chosen to provide the most 
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reliable extension of the source. As is apparent in a slice through the source, as shown in 
Figure [2] (bottom), a Gaussian profile might not be the correct representation of the source 
morphology, which seems to exhibit a tail extending mostly southwards. As can also be seen 
from the slice, the contamination of a possible point-source-like pulsar within the source is 
expected to cover one or two bins in the slice and therefore to contaminate the emission 
by ~ 20%. The source position determined using the XMM-Newton data (see Table [ID 
is compatible with both the ASCA source (at 16 h 40 s 42.24 s , -46°32'42.0" with a typical 
positional error of 1' on the ASCA position and a distance of 1.3' between the ASCA and 
the XMM-Newton best fit position) and the Swift XRT source (at 16 h 40 s 43.5 s , -46 31'38.6" 
with a positional error of 6" at a distance of 20"). As is apparent from Figure [21 the 
extended source XMMU J164045.4-463131 is positionally coincident with HESS J1640-465 
and located within the radio shell G338. 3-0.0, suggesting an association between the sources 
in these different energy bands. XMMU J164045. 4-463131 does not appear to be positionally 
coincident with the weak point-like radio source located within th e shell of G338. 3- 0.0 if the 
error on the radio position is indeed as small as several arcseconds (IGreen et al.lll999l ). Future 
radio studies are well motivated to identify the nature of this compact radio source. 



Table 1: X-ray sources detected in observations of HESS J1640-465 using the detection al- 
gorithm emldetect. The parameters given here are for the energy range between 0.5 keV and 
10 keV. The second column gives the name recommended by the XMM-Newton SOC and 
the IAU for source detections. Columns 3 and 4 give J2000 coordinates. Column 5 gives 
the error on the source position in arcseconds and column 6 gives the number of counts in 
EMOS1 and EMOS2 within a 35" integration region using events above 0.5 keV. 



Id 


XMMU J 


RA2000 


Dec 2 ooo 


Error RA 20 oo 


Counts 






(h:m:s) 


(d:':«) 


(") 




1 


164045.4-463131 


16:40:45.39 


-46:31:31.1 


8.5 


222 


2 


164029.6-462328 


16:40:29.58 


-46:23:28.0 


1.5 


153 


3 


164131.0-463048 


16:41:30.98 


-46:30:48.0 


0.6 


215 



A catalogue search did not yield any obvious counterparts to sources #2 and #3. 
Since the angular distance to these two sources is too large to be reasonably associated 
with the VHE 7-ray emission in HESS J1640-465, the following analysis will be focused on 
XMMU J164045. 4-463131. For the spectral analysis, the data was extracted using SAS 7.0. 
The reduced data set has been analysed and fitted using XSPEC (version 12.2.1). Several 
different background estimation techniques have been applied, including determining the 
background from the same field of view (both a ring around the source region and the source 
region mirrored at the centre of the camera) as well as from blank field observations. The 
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different background estimation techniques agree well within the statistical errors and there- 
fore in the following the background taken from a ring around the source in the same field 
of view will be used. The application of a background estimate from the same field of view 
has the advantage that all relevant backgrounds for this observation are included, such as 
the X-ray emission from the Galactic ridge, the residual particle background as well as the 
instrumental background. To encompass the whole extension of the source a source region of 
radius 75" has been used for the spectral analysis. Consistent results with larger error bars 
were achieved with source extraction regions of radius 50" and 100". Table [2] summarises the 
parameters of the spectral fitting. Two spectral models were fitted to the data, an absorbed 
power-law and an absorbed black-body spectrum. 



Table 2: XMM-Newton spectral properties of XMMU J164045. 4-463131 for a simultaneous fit 
to the EMOS1, EMOS2, and EPN data. The two fit models are an absorbed power-law and 
an absorbed black-body. The fit parameters for the absorbed power-law are the absorption 
density nn in units of cm -2 , the photon index T, and the normalisation at 1 keV (kik e v)- 
The fit parameters for the black-body spectrum are the absorption density nn again in units 
of cm -2 , the temperature kT in units of keV, and the normalisation at 1 keV (k^ev)- Also 
given is the integrated absorbed (i.e. observed) flux between 2 keV and 10 keV F 2 _io kcv (in 
units of erg cm -2 s" 1 for the two models. The errors given correspond to 90% confidence 
levels. 



Parameter 


Value (powerlaw) 


Value (black-body) 


EMOS1 Counts 




321 




EMOS2 Counts 




290 




EPN Counts 




807 




n H (cm" 2 ) 


(6.ii 2 i) x 10 22 




(3.6j&) x 10 22 


r 


1.741°;- 






kT (keV) 






1-563:3 


klkeV 


(2.6^) x 10- 4 




x 10- 5 


F 2 -io keV (erg cm" 2 s" 1 ) 


6.6 x 10~ 13 




6.4 x 10~ 13 


xVd-o.f. 


2702/3001 




2694/3001 



The power-law fit yields a hard photon index of 1.74 ± 0.1, the source is rather faint 
with a total detected (i.e. absorbed) flux in the 2-10 keV range of 6.6 x 10~ 13 erg cm -2 
s _1 . The column density nn = 6.1 x 10 22 cm -2 is marginally compatible with the average 
column density through t he Galaxy in this region (nc a iH = 2.2 x 10 22 cm -2 as reported by 
Dickey fc Lockmanl (119901 ) ) . The black-body model yields a lower column density tir = 
3.6 x 10 22 cm -2 . It should however be noted, that the fit parameters (tih and kT) are 
somewhat correlated. Both models fit the data well as shown in Figure [31 Comparing the 
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values of the reduced x 2 ; n o spectral model can be preferred over the other. While the results 
for the different background estimates are compatible within statistical uncertainties, the 
scatter might give a hint on the systematic error steming from the choise of the background 
and thus the systematic error was estimated to 0.2 on the photon index and 20% on the 
integrated flux. The result of the spectral analysis of the XMM-Newton data presented here 
is in agreement with the ASCA data, given the large error on the latter spectral results 
(photon index of 3.0+^, F 2 _ 10keV = 1.2 x 10~ 12 erg cm" 2 s~\ n R = 9.6tf;£ x 10 22 cm- 2 ). 
The spectral analysis shows however, that the photon index is rather hard compared to the 
ASCA spectrum. 



3. Interpretation of the Mult i- wavelength data 



The XMM-Newton data reveal the source XMMU J164045.4-463131 in positional coin- 
cidence with the unidentified ASCA source AX J1640. 7-4632. The source exhibits a rather 
hard power-law spectrum of photon index ~ 1.75, although a black-body spectrum of the 
X-ray emission with a temperature 1.56 keV, corresponding to 1.8 x 10 7 K provides an 
equally good fit. The source is relatively faint (-F2-10 kov — 6.6 x 10~ 13 erg cm~ 2 s _1 ) and 
shows a rather strong absorption of n H = 6.1 x 10 22 cm" 2 . XMMU J164045. 4-463131 is 
extended in nature with a compact core and a faint tail, resembling in morphology and spec- 
tral properties typical Pulsar Wind Nebulae (PWN). Asymmetric "trails" in PWN can be 
generated either a) through an asymmetric density distribution of the surrounding medium, 
preventing the expansion of the PWN on one side as e.g. seen in Vela X, b) dynamically by 
a supersonic motion of the pulsar with respect to the ISM, generating a bow-shock and a 
"cometary trail" (for a recent review of PWN see iGaensler fc Slandl2006l ). The hypothesis 
that XMMU J164045. 4-463131 is a PWN (stemming from the hard energy spectrum and 
the morphology) is strengthened further by the location of this object within the bound- 
aries of the shell-type radio SNR G338. 3-0.0. This PWN scenario could finally be confirmed 
by the detection of a pulsar coincident with XMMU J164045. 4-463131. No pulsation has 
been found in the current XMM-Newton data and no know n radio pulsar is located towards 
XMMU J164045. 4-463131 in the current ATNF catalogue jManchester et all 12005 1). Future 
deep radio or X-ray timing observations might serve to detect the associated pulsar. The 
"beaming fraction", i.e. the faction of 4 ir steradians covered by the pulsar beam during 
one rotation, corresponds to the probability that the beam sweeps the line-of-sight of an 
observer. So far, no agreement has been reached on the relation between the beaming frac- 
tion and the period of the pulsar, and estimates range from beaming fractions of 30% to 



100% for a typical 100 ms pulsar (INarayan &: Vivekanandl Il983l ; iLyne &: Manchester! 11988 



Tauris &: Manchester! Il998l ). In any case, the potential pulsar to the PWN discussed here 
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might be beamed away from our line-of-sight. In this case, the confirmation of the PWN 
picture must come from the detection of electron cooling, resulting in a softening of the X-ray 



spect rum away from the pulsar as seen in many other PWN (e.g. in G 21.5-0.9. ISlane et al. 



2000|)). Given the faintness of XMMU J164045.4-463131 the current data do not allow for 
the detection of such an effect. Very deep high angular-resolution Chandra studies will finally 
confirm or rule out this scenario. Nevertheless, the XMM-Newton data are very suggestive 
of a PWN, since the hard X-ray spectrum, the extended nature, the position within an radio 
SNR and the spatial coincidence with a VHE 7-ray source seem unlikely to occur by chance. 
In the following discussion we will therefore assume that the emission in the radio, X-rays 
and VHE 7-rays are connected and characterise a composite SNR. 



(Aharonian et al. 


2005a. 


2006a: 


Funk et al. 


2006) 



coincidence of an extended TeV 7-ray source with an extended hard spectrum X-ray source 
located within the radio shell of an SNR. As in the case of HESS J1813— 178 it is not clear, 
whether the 7-ray emission originates in the shell of the SNR, or from the core, since the 
angular resolution of VHE 7-ray instruments is too coarse to resolve the small structures 
seen in X-ray and radio observations. The "Gaussian equivalent width" of the radio SNR 
can be defined as o"snr = (o~l + Oy — Smooth) ^ 2 ( a x an d cr y the Gaussian width in RA and 
Dec direction of the radio emission region smoothed with a Gaussian of size r smoot h = 1.2' as 
used in the H.E.S.S. analysis). This "equivalent width" of the radio shell of G338. 3-0.0 is 
2.5' ± 0.2', perfectly compatible with the size of the VHE 7-ray emitting region with an rms 
of 2.7' ± 0.5'. As in the case of HESS J1813-178 both the radio shell and the central X-ray 
core must therefore be considered as viable VHE 7-ray emitters. However, there are several 
notable differences to HESS J1813-178: 

• The X-ray emission of HESS J1640-465 is much weaker. This has several important 
consequences: any hard X-ray counterpart falls well below the sensitivity limit of 
current detectors such as INTEGRAL and Suzaku. Thus no hard X-ray counterpart to 
XMMUJ164045.4-463131 exists as seen for HESS J1813-178 (suggesting acceleration 
to PeV energies). The faint X-ray emission also does not allow for an exclusion of 
a thermal origin of the keV emission. Future deep X-ray observations might finally 
distinguish between a power-law or a black-body emission model. 



An unidentified EGRET source (3EG J1639-4702) flHartman et al.l Il999f ) is in posi- 
tional coincidence with G338. 3-0.0. If the sources are associated, this increases the 
7-ray power even further. 

So far, no dedicated search for a radio pulsar has been performed within the remnant 
and thus so far no pulsar has been found in G338. 3-0.0. As discussed above, the 
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pulsar light-cone might not sweep us as observers. There is however a weak radio 
point-source located at 16 h 40 m 47.75 s , — 46 d 32'03.3" within the remnant, that could 
possibly be identified as a pulsar or even a "relic" PWN in deep observations. This 
source, which can be seen in Figure [H lies ~ 40" from XMMU J164045. 4-463131. 

• The distance to HESS J1640-465 (~ 8.6 kpc) is only poorly known, constrained only 
by the E — D relation for G338. 3-0.0 (with E being the surface brig htness and D 
the distance). Such estimates are notoriously unreliable (IGreenl 120041 ) . The nearby 
HII- regi ons shown in Figure [1] hav e measured radial velocities and implied distances of 
3-4 kpc (IWhiteoak fc Greenlll996l ). and G338. 3-0.0 may be associated with this region 
of star formation. The relatively strong H.E.S.S. and EGRET fluxes (provided that 
these objects are associated) would also argue for a distance <C 10 kpc on energetics 
grounds. Nevertheless, in the following a distance of ~ 8 kpc will be assumed. 



Whilst several morphological properties of HESS J1640-465 are very similar to HESS J1813- 
178, spectral properties differ quite strongly and the first two differences listed above ar e 
reminiscent of another VHE 7-ray source: HESS J1825-137 (jAharonian et al.ll2005a . 12006a ). 
This object is the only kno wn VHE 7-ray source to exhibit energy dependent morphol- 
ogy JAharonian et al.ll2006ch . For HESS J1825-137, an X-r ay PWN was found in XMM- 
Newton observations of the energetic pulsar PSR J1826-1334 (IGaensler et al.ll2003l ). extend- 
ing asymmetrically to the south of the pulsar. Similarly, the VHE 7-ray emission shows an 
asymmetric emission, extending to the south of the pulsar, however on a completely different 
scale than the X-ray emission (the X-ray emission extends 5", whereas the 7-ray emission 
extends ~ 1° to the south). Taking into account the different loss timescales of the 7-ray 
and X-ray emitting electrons, the situation can be plausibly explained by the 7-rays being 
generated by lower energy electrons than the X-rays. This scenario was strengthened by the 
detection of energy dependent morphology in the VHE 7-ray sour ce, leading to a convinc - 
ing identification of this object as the PWN of PSR J1826-1334 JAharonian et alJl2006ch . 
PSR J1826-1334 is a powerful "Vela-like", ~20 kyr old pulsar. For such evolved pulsars, 
in which the spin-down luminosity has changed substantially during the lifetime, braking 
effects are no longer negligible and a decrease of the parti cle injection rate with age must 
be taken into account. According to lPacini fc Salvatil (119731 ). the spin-down luminosity E(t) 
evolves in time as 



E(t) = E(0) 



n+l 



with to 



Po 

Pn(n - r 



(with n the braking index and To the initial spin-down timescale of the pulsar). If HESS J1640- 
465 indeed represents a PWN similar to HESS J1825-137, a relatively powerful, 10 4 — 10 5 
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year old pulsar should exist within G338. 3-0.0. The age of G338. 3-0.0 has been estimated 
to ~ 2 x 10 4 years, assuming expansion in the Sedov phase into a medium of density 10 
cm -3 , and using the distance estimate of ~ 8 kpc. 

Taking into account the examples of both HESS J1813— 178 and HESS J1825-137, in 
the following the multi-frequency properties of HESS J1640-465 will be discussed with the 
aim to connect the radio, X-ray and VHE 7-ray emission into a consistent picture. Table [3] 
summarises the properties of the X-ray and of the 7-ray emission. 



Table 3: Comparison between size and spectral properties of XMMU J164045. 4-463131 and 
HESS J1640-465 for an assumed distance of 8 kpc. 





XMMU J164045.4-463131 


HESS J1640-465 


Angular Size 


0.45' ± 0.05' 


2.7' ±0.5' 


Linear Size (pc) 


~ 1^8kpc 


~ 6d 8kpc 


Apparent Luminosity (xlO 33 erg/s) 


du 8kpc 




Photon Index 


1.75 ± 0.12 stat 


2.42±0.15 stat 



Figure [5] compares the spectral energy distributions of HESS J1640-465 and HESS J1825- 
137. The main characteristics that have to be explained by the modelling are that the VHE 
7-ray source is an order of magnitude larger than the X-ray source and that the VHE 7-ray 
power is more than an order of magnitude larger than the X-ray power. As in the case of 
HESS J1825-13 7, it is rather difficult to estab lish a multi- wavelength connection in a leptonic 
one-zone model (lAharonian fc Atoyan 1119991 ) with constant injection and a single population 
of electrons responsible for the whole emission. For plausible radiation and magnetic field 
energy densities, such models generally predict X-ray synchrotron fluxes equal to or greater 
than the associated TeV IC fluxes, unless suppress ed magnetic fields or rathe r peculiar en- 
vironments like in the Galactic Centre are invoked (IHinton fc Aharonianll2006l ). As it seems 
hardly possible to find satisfactory parameters for such a simple model, a time- dependent 
rate of injection of electrons into the nebula will be considered. For the following argument, 
the existence of a ~ 2 x 10 4 year old pulsar produced in the supernova explosion of G338.3- 
0.0 will be assumed. For younger ages, no satisfactory fit to the data can be found in the 
frame of the model described here. It will be further be assumed that the size of the X-ray 
source is limited by electron synchrotron cooling which in turn limits the age of particles seen 
at X-ray energies. To explain the different sizes of the emission regions in X-ray and VHE 
7-rays, cooling times of the X-ray emitting electrons must be short in comparison to the age 
of the source, while cooling times for the VHE 7-ray emitting electrons must be long in com- 
parison to the age of the source. This situation is rather likely for a system of this age: For a 
typical 10 /xG field, 4 keV synchrotron X-rays are generated by ~ 100 TeV electrons, which 
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lose their energy on timescales of ~1,200 years. In contrast, 0.8 TeV IC 7-ray are emitted 
by ~ 3 TeV electrons, which cool on timescales of ~ 32, 000 years. A natural explanation 
for the different sizes of the emission regions is to invoke different populations of electrons 
with different cooling timescales and hence propagation distances. In this scenario the size 
of the X-ray nebula is limited by electron cooling on ~ 1000 year timescales, in contrast 
the angular size of the 7-ray source reflects the propagation speed of electrons injected soon 
after the birth of the pulsar ~ 2 x 10 4 years ago. 

The curves in Figure [5] illustrate the model scenario invoked to connect XMMU J164045.4- 
463131 to HESS J1640-465 for populations of electrons injected at different slices in time t 
during the lifetime of the pulsar. A canonical E~ 2 injection spectrum was used with a lower 
energy cut-off at 10 GeV and an exponential cut-off at 1 PeV, as the rather hard X-ray spec- 
trum yields a lower limit for the cutoff energy of a few hundreds of TeV. A time-independent 
magnetic field strength of 10 fiG was assumed for the whole nebula. The injection rates of 
electrons varies with time, assuming that the power injected into electrons follows the spin- 
down power of the pulsar which varies according to Equation [1] (with n = 3 and To = 300 
years). The dotted line shows the electrons injected in the first 2000 years after the birth of 
the pulsar (20,000 years ago), i.e. "old" electrons. The injected electron population evolves 
in time, taking into account synchrotron and IC energy losses. The population of "old" 
electrons totally dominates the synchrotron emission below X-ray energies and is responsible 
for the VHE 7-rays detected in HESS J1640-465. The populations of "youngest" electrons 
(dashed curve) injected in the last 2000 years are less drastically cooled, as can be seen in 
FigureO According to the model shown here, these youngest electrons are responsible for the 
bulk of the X-ray emission as seen in the compact X-ray source. To explain the X-ray source 
as synchrotron emission from a PWN, the required (present day) luminosity in accelerated 
electrons is 1.4 x 10 36 (B/10yuG)~ 2 (d/8kpc) 2 erg/s, or about 50% of the spin-down power of 
a pulsar like PSR J1826-1334 (3 x 10 36 erg/s), assumed to be a representative of the pulsars 
powering these VHE 7-ray PWN systems. The total electron luminosity integrated over the 
lifetime of the pulsar amounts to 1.4 x 10 48 ergs. An upper limit from the 8 /zm GLIMPSE 
Spitzer data, integrating within the area of the VHE 7-ray source, yields an unconstraining 
upper limit of ~ 1 x 10~ 8 erg cm~ 2 s _1 (at ~ 0.2 eV). It should however be noted that in the 
derivation of this limit no subtraction of prominent IR sources has been applied. Although 
the figure shows how to accommodate the different sizes and fluxes of the X-ray and VHE 
7-ray production regions, it must still be seen as an oversimplified picture. A realistic model 
requires the inclusion of the time (and space) dependence of the magnetic field, as well as a 
proper treatment of the propagation of particles within the nebula. 

Also apparent from these model curves is that it is rather difficult to connect HESS J1640- 
465 and XMMU J164045. 4-463131 with the unidentified EGRET source 3EG J1639-4702 
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within the framework of an IC scenario. However, an association of 3EG J1639-4702 and 
HESS J1640-465 is rather natural in two alternative scenarios: a model in which the VHE 
7-rays are generated by electron Bremsstrahlung or in a hadronic emission model. If the ISM 
density in this region is high (n is more than about 100 cm -3 ), then Bremsstrahlung may be 
the dominant 7-ray emission mechanism at TeV energies for a population of ultra-relativistic 
electrons. If HESS J1640-465 and 3EG J1639-4702 are indeed associated, Bremsstrahlung 
may be the most natural mechanism to explain the high energy part of the SED. 

Alternatively, the EGRET and H.E.S.S. detected 7-rays could be generated by hadronic 
interactions of accelerated protons via 7r°-decays. In this case the 2-10 keV X-ray flux can 
be connected to these accelerated hadrons by secondary electrons generated in the hadronic 
interactions. Assuming a distance of 8 kpc, the measured 7-ray flux above 200 GeV sug- 
gests a total energy in the accelerated hadrons of W p = 7.6 x 10 50 <i/(8 kpc)n/(lcm~ 3 ) 
erg. This energy represents a plausible injection from a single SNR for values of n more 
than 1 cm -3 and/or distances significantly closer than 8 kpc. However, in this scenario, 
secondary electrons would be produced in the same region as the TeV 7-rays, and similar 
sizes for the X-ray and 7-ray sources might naively be expected. The observed factor of six 
difference in size therefore counts against an association of XMMU J164045. 4-463131 and 
HESS J1640-465 in this scenario. Whilst a hadronic scenario for 7-ray emission remains 
viable, if XMMU J164045. 4-463131 and HESS J1640-465 are unrelated, then the excellent 
positional agreement between these two sources can be considered somewhat surprising. 

Given the fact that there are several other prominent candidates for 3EG J1639-4702 
(as described above) and that the EGRET source is the least well measured object in this 
context, it seems that a strong claim of a possible connection has to await the launch of the 
upcoming GLAST satellite. The dashed line in Figure [5] shows the 1-year sensitivity of the 
GLAST-LAT, including the diffuse Galactic background and the instrumental background. 
The curve demonstrates that the LAT is well suited to shed new light on this region in the 
MeV to GeV range via its superior angular resolution (and sensitivity) in comparison to 
EGRET. 



4. Summary and conclusion 

The detailed XMM-Newton X-ray data taken towards HESS J1640-465 show a hard- 
spectrum extended object towards the centre of the SNR G338. 3-0.0, suggesting the detec- 
tion of a new composite SNR in which the X-rays are generated by synchrotron emission from 
a central PWN. The ultimate proof of this scenario requires the detection of a coincident 
pulsar, in either X-ray or radio observations, or the detection of spectral cooling character- 
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istic of PWN systems. Both the shell of the SNR and a central PWN are viable VHE 7-ray 
emitters and the two scenarios cannot be distinguished at this point (the situation is similar 
to the case of HESS J1813-178). If the VHE 7-ray are generated by inverse Compton emis- 
sion associated with a PWN, the spectral energy distribution shows an interesting similarity 
to another VHE 7-ray PWN HESS J1825-137. As is the case for that object the low ratio of 
the X-ray power to the VHE 7-ray power in HESS J1640-465 suggests a time dependent rate 
of injection for the relativistic electrons responsible for the X-ray emission and an older (and 
more numerous) population of electrons producing the VHE 7-ray emission. IR to soft X-ray 
measurements of the synchrotron nebula associated with HESS J1640-465 could confirm this 
picture. However, any synchrotron emission from HESS J1640-465 in this energy range ap- 
pears to be buried beneath thermal emission. An extension of the 7-ray spectrum to higher 
energies (for example from a deeper observation with H.E.S.S.) would allow us to probe the 
IC emission of the electrons responsible for the > 1 keV synchrotron emission. The detection 
of the extended X-ray source in the centre of the Supernova remnant G338. 3-0.0 following 
the detection of a VHE 7-ray source demonstrates that X-ray follow-up observations of TeV 
7-ray are well motivated and provide important insights into the nature of these objects. 
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Fig. 1. — The XMM-Newton field of view for the observation on HESS J1640-465. This 
composite image of the MOS-1 and MOS-2 count maps shows the energy range below 2 keV 
(top) and above 2 keV (bottom), the counts were adaptively smoothed. Three prominent 
X-ray sources can be seen lining up with previously reported ASCA and Swift XRT sources 
(t he cyan circ l es den ote the position of the Swift sources #1, #2, and #3 as reported 



in 



Landi et al.1 (120061 )). Source #1 is coincident with the known X-ray source AX J1640.7- 
4632 located within the radio shell of G338.3-0.0. The 843 MHz MOST radio data are 
shown as white contours. The weak central radio object within the shell is also apparent in 
these contours. Also visible from this image is stray-light contamination from the close-by 
low-mass X-ray binary GX 340+0. 
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Fig. 2. — Top: Zoomed view on the XMM-Newton field of view for the observation on 
HESS J1640-465 for events above 2 keV (composite image of the MOS-1 and MOS-2, adap- 
tively smoothed). The white contours denote the 843 MHz MOST radio data showing the 
inner edge of the radio shell of G338. 3-0.0. The black square along with the solid cyan 
lines denote the best fit position (and la error) of the VHE 7-ray source HESS J1640-465, 
the dashed cyan circle indicates the rms extension of HESS J 1640-465. Also shown are the 
XMM-Newton best fit position as detected with emldetect as well as the best fit position of 
the Swift XRT X-ray source. Bottom: Slice along the dotted red box (top Figure) through 
the X-ray (black) and radio (blue) emission. The difference between the shell-like structure 
in the 843 MHz radio data and the compact core with extended tail towards the south in 
the XMM-Newton data is clearly apparent. 
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Fig. 3.— XMM-Newton X-ray energy spectrum of XMMU J164045.4-463131 for two different 
spectral models, an absorbed powerlaw (top) and an absorbed black-body (bottom). All 
detectors have been simultaneously fit (EMOS1: black, EMOS2: red, EPN: green) for an 
extraction radius of 75". The plot shows flux points with a minimum significance of 5a. 
The straight lines shows the combined fit to the data, folded with the instrument response 
function of the corresponding detector). The lower panel of each plot shows the residuals of 
the fit, illustrating the good match of both fit functions to the data. 
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Fig. 4. — Spitzer Space telescope data of the GLIMPSE survey at 8 fim ([Benjamin et al. 



in 



200 31) along with 843 MHz MOST radio data shown as white contours. As described 



Whiteoak fc Greed (119961 ). the so-called radio-bridge extending to the north of G338. 3-0.0 
coincides very well with a bright Hll-region in which dust emission shows ongoing star for- 
mation. This region is located at a distance of ~ 3 kpc and might have played an important 
role in the development of the SNR if both objects are located at the same distance. Also 
shown are the position and extension of HESS J1640-465 (pink) and the combined MOS1 
and MOS2 XMM-NewtonX-ray contours (green) all coinciding with the radio SNR G338.3- 
0.0 described in the text. Also shown are the known Masers in the region, indicating an 
interaction between shock waves and the ambient medium and illustrating that the region 
is still active in terms of star formation. 
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Fig. 5. — Spectral energy distribution of HESS J1640-465 for the model with a time de- 
pendent rate of injection of relativistic electrons as described in the text. For comparison 
the spectral data for the VHE 7-ray PWN HESS J1825-178 are shown as grey bands. The 
XMM-Newton data between 2 and 10 keV are corrected for the absorption and a systematic 
error band of 0.2 on the spectral index and 20% on the flux level has been added. The 
different curves indicate emission from electron populations injected at different time slices 
in the age of the pulsar (assumed to be 20,000 years of age). The injection rate was assumed 
to vary in accordance with the change in spin-down power as described in the text. Youngest 
electrons (injected in the last 2000 years) responsible for the X-ray emission are shown as 
dashed lines, oldest electrons (injected in the first 2000 years after the pulsar's birth) are 
shown as dotted lines, dominating the VHE 7-ray emission. The red curves show the Inverse 
Compton component, whereas the blue lines show the synchrotron emission in the chosen 
model. The radiation fields for the inverse Compton emission are assumed to be nominal 
Galactic radiation fields taken from lPorter et al.l (120061 ). The model parameters are specified 
in the text. The black dashed line shows the 1-year sensitivity curve for the GLAST-LAT 
taking into account the diffuse Galactic and the residual instrumental background. 



